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We study a doped antiferromagnet (AF) using a rotating reference- frame. Whereas in the labo- 
ratory reference-frame with a globally fixed spin-quantization axis (SQA) the long-wavelength, low- 
energy physics is given by the 0(3) non-linear cr-model with current-current interactions between the 
fcrmionic degrees of freedom and the order-parameter field for the spin-background, an alternative 
description in form of an U(l) gauge theory can be derived by choosing the SQA defined by the 
'• local direction of the order-parameter field via a SU(2) rotation of the fermionic spinor. Within a 

■ large-iV expansion of this U(l) gauge theory we obtain the phase diagram for the doped AF and 

| identify the relevant terms due to doping that lead to a quantum phase transition at T = from 

the antiferromagnetically ordered Neel phase to the quantum-disordered (QD) spin-liquid phase. 
Furthermore, we calculate the propagator of the corresponding U(l) gauge field, which mediates 
a long-range transverse interaction between the bosonic and fermionic fields. It is found that the 
strength of the propagator is proportional to the gap of the spin-excitations. Therefore, we expect 
as a consequence of this long-range interaction the formation of bound states when the spin-gap 
£f*) , opens, i.e. in the QD spin-liquid phase. The possible bound states are spin-waves with a (spin-) gap 

1 in the excitation spectrum, spinless fermions and pairs of fermions. Thus, an alternative picture for 

charge-spin separation emerges, with composite charge-separated excitations. Moreover, the present 
treatment shows an intimate connection between the opening of the spin-gap and charge-spin sepa- 
ration as well as pairing. 
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lO ' I. INTRODUCTION 

On ' 

Recently much attention has been focused on the problem of understanding the magnetic as well as transport 
properties of high-temperature superconductors (HTSC) as a function of doping. The undoped parent compounds like 
La2CuC>4 are antiferromagnetic insulators with 5=1/2 local moments on the copper-sites Upon doping charge- 
carriers into the CuC>2 planes, the long-range magnetic order is destroyed and these materials become superconducting 
at low temperatures, while in the normal state unconventional properties are observed . The particular interesting 
question arising here is whether these anomalous features are directly related to the proximity to the magnetic quantum 
phase-transition or not. 

Since in the relevant low-doping regime the spin-spin correlation length £ is large Q , it is appropriate to address this 
question from a field-theoretic point of view. Therefore we will study this problem within an effective long- wavelength, 
low-frequency theory, which can be derived in a rigorous way from a realistic microscopic model, i. e. without semi- 
classical approximations and phenomcnological assumptions. By means of a SU(2) rotation of the fermionic spinor 
one can further derive two alternative descriptions of the doped AF [Q. In the laboratory reference-frame with a 
globally fixed SQA, the doped AF is described by the 0(3) non-linear cr-model for the order-parameter field n of the 
spin-background and current-current interactions between the doped fermions and the order-parameter field. On the 
other hand, in the rotating reference-frame with local SQA defined by the local direction of n, the long-wavelength 
low-energy theory for the doped AF is given by an U(l) gauge theory with bosonic and fcrmionic degrees of freedom 
minimally coupled to an U(l) vector gauge field. The present work presents an analysis of the gauge theory in the 
rotating reference- frame. 

Within a large- TV expansion of the U(l) gauge theory we will explicitly show 1) that there occurs a quantum phase- 
transition from the antiferromagnetically ordered Neel-state to a QD spin-liquid phase at T = induced by doping 
and 2) that the U(l) gauge field mediates a long-range interaction between the fermionic and bosonic fields, where the 
strength of the coupling is given by the spin-gap. Therefore, bound states between excitations with opposite charge 
(as defined by the coupling to the gauge field) are expected to form in the QD phase. The possible bound states 
correspond to spin- waves with a gap in the excitation spectrum, spinless fermions and pairs of fermions. In this frame 
the phenomena of charge-spin separation and pairing are intimately connected with the opening of a spin-gap and 
have the same physical origin, i. e. the long-range interaction induced by the vector gauge field. 

In order to clarify the connection of the gauge theory studied here with a microscopic model, we summarize in the 
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following the main steps involved to reach the effective field-theory Q . Our starting point is the strong coupling limit 
of the three-band Hubbard model for the Cu02 planes, which is given by the spin-fermion model 0] and acts in the 
subspace of singly occupied Cu-sites: 



H 



E 



Jk^^Ri ' Si + Jh /J Si ■ Si 



(1) 



and Cj t u create and annihilate holes with spin-projection a =|,| in an oxygen-orbital. The kinetic term 



where c] CT a ^ oj i<r 

describes hopping processes between oxygen-sites and can include a direct oxygen-oxygen hopping as well as an effective 
hopping via the central copper-site. Furthermore the Hamiltonian contains a non-local, Kondo-like interaction between 
the localized Cu-spins Si and the neighboring O-holes Ri — Ylti.i*-.*) (— l) aij+aik c) cr a p Ck,p (here the summation 

a,0 J ' 

goes only over the indices j and k denoting the oxygen-sites around the central copper-site i), where u are the Pauli 
matrices, as well as an antiferromagnetic Heisenberg superexchange interaction between nearest neighbor Cu-spins. 
For a connection of the parameters of the spin-fermion model with the ones of the three-band Hubbard model as well 
as the definitions of the phase factors ay which accounts for the d- and p-wave symmetry of the underlying orbitals 
we refer to Ref. [Q. 

In order to construct an effective theory for the long-wavelength low-energy sector of the spin-fermion model we 
have to integrate out all the high-energy modes of the model. For the spin part this is done by means of an expansion 
around the adiabatic limit within the method of generalized Berry phases ^|||. The continuum-limit in the order- 
parameter field (fii of the spin-background is achieved by means of a gradient-expansion around the AF saddle-point 
<p i = (— where the lattice constant a (or the inverse of an equivalent ultraviolet cutoff) serves as the small 

expansion parameter. Consequently, the Kondo-like interaction in the Hamiltonian Eq. ([!]) goes over to essentially 
the same form but with the order-parameter rii instead of the spin-operator Si. Therefore it becomes a staggered 
potential for the fermions in the antiferromagnetic ground-state and together with the hopping parts it determines 
the dispersion of the fermions which has minima located at k m i n = (±7r/2, ±7r/2). Since we are interested in the low 
doping regime and low-energy sector of the model it is justified to substitute the whole band structure by parabolic 
bands located at k min . In this way we obtain the following continuum theory: 





S N i a = ^- J dr I <Hr 
o 



[ P dr [ d 2 r Sd T c + -L d r J d r c - lr (d r n) 2 Sc + V 7(1 ^-jJ 



S F = dr d 



H = (r, t) 



(2) 
(3) 



where the first two terms in Sf describe the kinetics of the fermions with perpendicular m± and parallel my masses 
with respect to the borders of the magnetic Brillouin-zone. The remaining parts denote the interaction of the fermions 
with the antiferromagnetic spin-background, whose dynamics is controlled by the 0(3) non-linear cr-model with 
g = 2\[2 a/s being the coupling constant and c = 2\/2 aJns being the spin-wave velocity, s denotes here the length 
of the spin. The third term of Sf gives a renormalization of the spin-stiffness of the non-linear cr-model, whenever the 
corresponding site is occupied by a hole. The fourth term gives current-current couplings similar to the ones obtained 
by Shraiman and Siggia for the t-J model 0. For the temporal component we have a coupling of the local spin- 
density = c^crc of the fermions to the background magnetization = in x d T n whereas the spatial components 
denotes a coupling between the spin-current of the holes = i{d r c^crc — c^crd r c) and the magnetization current 
of the background = n x d r n. The connections of the corresponding coupling constants as well as the masses 
to the microscopic parameters of the spin-fermion model are given in Ref. 0]. For simplicity we will use isotropic 
parameters 7 = 7± = 711 and m = mi = TO|| in the following. However, such a situation is close to the case given by 
a parameter set for the underlying three-band Hubbard model, which, as quantum Monte Carlo simulations showed 
M, is appropriate for a quantitative description of various physical quantities of the HTSC. 

Alternatively one can describe the system using a reference-frame, which is defined by a local SQA for the fermionic 
spinor in the direction of the order-parameter field n. This change of reference-frame is defined by the following local 
SU(2) rotation of the fermionic spinor 



p(r, t) — U\r, t) c{r, r) 



(4) 



where the SU(2) rotation U must fulfill U^r^r 
rotation matrix U 



cr ■ n(r,r) U(r,r) = a z . Using a CP 1 representation [B| for the 
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Zl -Z 2 
Z2 Zl 



(5) 



with ZZ = 1 and Z = (z\, z<i) denoting a bosonic spinor, the transition to the rotating reference-frame is achieved by 
the following set of equations: 

uniform SQA < — > rotating SQA 

< > Pt* 

d„ + iA^ z (6) 



<9 M - -cr • (n x <9 M n) 
-cr • (n x <9 M n) 



Here we have introduced a composite gauge field and denoted the off-diagonal contributions by B^ and 



An = -iZd^Z 



B„ 



i {zid^z 2 - z 2 d^zi) = —Za v d,j,Z T , 



(7) 
(8) 
(9) 



which are related to the SU(2) rotation via Wd^U — icr z A fl + iB^<j + + iB^a . As a result of this transformation we 
get the following fermionic action 

S F = dT j ST \f ( D r + 2 7r A' T ) p + ±-D F r pt D F p + 7 [2 (d r p^K r p - p^Krdrp) - ^ K r K rP ] J , (10) 

with D F — dft+ia 2 A^, (/j G (r, r)) denoting a covariant derivative for the fermions and 7 r = 7r + 1 and 7 = 7+ l/2m. 
Accordingly the 0(3) non-linear cr- model transforms within this scheme by applying the Hopf projection pi defined 
by 



n = ZaZ , n 2 = 1 ZZ = 1 . 



to Eq. (§) into the CP 1 model 

Sep* = 



9 Jo J 



d z r <j DfZ D^Z + \r>®Z Z 



(11) 



(12) 



where = <9 M — iA^ denotes a covariant derivative for the bosons. Here and in the rest of the paper the superscripts 
F and B refer to the contributions of the bosons and the fermions, respectively. The whole action in the rotating 
reference-frame S = Sf + Scp 1 is invariant under a t/(l) gauge transformation. 

From the definitions of covariant derivatives we get the following charges of the bosons and fermions with respect 
to the gauge field A^: 





charge 


Zl 


+ 1 


Z2 


+ 1 


P\ 


+ 1 


Pi 


-1 



(13) 



The gauge theory defined by Eqs. (|T(]) and (jlj) is the starting point for the considerations in the present work. 

The outline of the paper is as follows. Section [H] is devoted to a large- A expansion of the U(l) gauge theory, where 
we will show 1) that the theory describes a doping induced qua ntum phase transition from an antiferromagnetically 
ordered state to a quantum disordered spin-liquid state (Sec. [I A ) and 2) that the transverse gauge fields that mediate 
a long-range inte racti on between the bosons and fermions has a strength given by the mass of the Z-bosons, i.e. 
the spin-gap (Sec. [I B ). In Section [II we will discuss the physical consequences of the gauge field and consider the 
connection between the opening of a spin-gap and charge-spin separation and pairing. In this section we will also 
address the implications of the bound states on various physical properties of the cuprates in the low-doping regime 
and we will give a brief summary of our main results and conclusions. A short account of the results above is given 
in Rcf. HI. 
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II. LARGE- TV EXPANSION FOR THE GAUGE THEORY 



The gauge field A M is not an independent dynamical degree of freedom since it is composed by two Z-bosons 
and has no kinetic part in the action (|l0|). Nevertheless, as usual such a kinetic term will arise due to the 
fluctuations of the bosonic as well as of the fermionic components. The calculation of such fluctuations is, however, 
not straightforward since the coupling constant for the gauge fields is of order unity, implying that there exists no 
small expansion parameter in the theory. In order to overcome these difficulties a large- N expansion is implemented, 
where N is the number of bosonic or fermionic flavors [ jl2| . 

Within the large- N expansion it will be shown that doping drives the system from an antiferromagnetically ordered 
state at T = to a quantum disordered spin-liquid state. A mass for the Z-bosons is generated dynamically, such that 
a " spin-gap" enters naturally in the description of the low-energy properties of the system, much in the same way as 
in recent analysis of the 0(3) non-linear cr-model In the present case, however, the non-diagonal contributions 

oc will lead to a doping dependence of the mass and, hence, a doping induced quantum phase-transition takes 
place. In the course of the large- TV expansion also the dynamics of the gauge fields is obtained, where the new energy- 
scale generated dynamically (i.e. the mass for the Z-bosons) determines the strength of the propagator of the gauge 
fields. This fact reveals an intimate connection between the spin-gap and the properties of the charge-carriers in the 
QD phase. 

As a first step we introduce L copies of fermions with spin ] and with spin J, 



(14) 



Accordingly the CP 1 fields are generalized to CP N 1 ones 



(zi,z 2 ) — ► Z = (zi, ...,Zn) 



(15) 



with TV = 2L so that the ratio of fermionic and bosonic degrees of freedom remains the same after this generalization. 

We proceed further by introducing sources for the fields and obtain (after proper renormalization of the fields) the 
following generating functional 



z[j, J,f},fiM = J vzvzvpVp^ J] s (zz- 2 -^ 



■ exp 



s, 



CP" 



J dr J d 2 r (fjp + p*n + ZJ + JZ+ T^QuA 



(16) 



Since the action Sf is bilinear in the fermionic variables they can easily be integrated out: 



S F = -iVTrln A F = -WTrln(A^ + Ap + Ap) 



(17) 



Here we have divided the fermionic contribution into three terms 
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A A 

^F 



2m r 



2gd T 
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2m 



ia z 2—d r 

0Q r 



+ 72(^) [2 (B^ 



ia z d r 
YB r a A 



S 



2g ia z 



1 S 2 



SQ r 2m 8Q 2 r 



+ 7(^) *B r B r 
d r + d r (B r <j- + B r a A 



(18) 
(19) 

(20) 



the first one is the free kinetic part, the second one contains the contributions from the gauge field A^, which we 
have substituted by a functional derivative with respect to its source ((g/2N) A^ — > as is clear from Eq. JT6|), 

and the third term contains the off-diagonal contributions K^. This division turns out to be useful, since due to 
the underlying SU(2) structure the A^ and -K^-terms in Sf do not mix within the large- N expansion and we can 
therefore treat them independently form each other. 

We first show that the K^-terms modify the CP 1 ^^ 1 model. An expansion of the non-diagonal part in powers of 
l/VN yields: 
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K 



d 3 pN 



(—)' 

\2NJ 

(J-Y 



AB T (jp)B s {-p) 



+ I d3 P N {^) *B T (p)B T (-p) 2 



d\G F {q) + -f / d 3 qG F (q)(2q+p) r (2q + p) s G F (q+p) 



d 3 qG F (q)G F (q+p) , 



(21) 



where G F (q) — \/(iv n — q 2 /2m) is the fermionic Greens function and where we have introduced the short hand 
notation J d 3 q = I d 2 q/(2n) 2 . Evaluating the integrals in Eq. ( |2l| ) in the limit p — > and u/\p\ — > (after 

proper analytic continuation) and noticing that (g/2N)B fJi B u = D^Z Z one gets for S F together with the pure 
bosonic contribution the following generalized CP 1 ^ -1 model: 



gen. _ „ qK _ 

(JpN-1 — OCPN-I ■+■ Dp — 



dr I d 2 r I ? r D*Z Df Z + —f S T D^Z Z 



with doping dependent coupling constants 

/* = 1 + 2g S (7 - 2m 7 2 ) = 1 - 2gS (1 + 2 7 to) 7 
^ = l + 2 5 c 2 7^e( ( 5), 



(22) 

(23) 
(24) 



where 0(<5) = 1 for S > and zero otherwise. Before we can integrate over the bosonic variables in the path integral, we 
have to decouple the quartic terms. This is achieved by means of the following Hubbard-Stratonovich transformation 



eXP {^i dT J ^[-(ZBMZd^+iQuiZdpZ)] 
X>A M exp 




dr I d 2 r 



2 2i 1 

— j-A^A^ H — j=X^(Zd^Z) H — ^=A M Q M 



iV 



5/^ 
87V 



QuQl 



(25) 



where we have rescaled the bosonic fields in such a way that now the constraint reads S (^ZZ — -^re- 
absorbed the factor ff/c 2 in a redefinition of the inverse temperature (c/y/*)/3 — > /3. Since the introduced Hubbard- 
Stratonovich field A M couples now linearly to the source it plays from now on the role of the gauge field. Furthermore 
we include the constraint for the bosonic variables via a lagrange multiplier field a into the action: 



and we have 



27V 



T>a exp 







dr / d 2 r 



N 



2N 
Jltt 



(26) 



After these modifications the action is bilinear in the bosonic variables and they can be integrated out. In this way 
we get the following effective action 



S cS = NTr In A B - NTr In (A% + A F ) + f dr f 

91 Ir JO J 



d 2 r a{r 1 t) 



with 



A« = -D°D 



B n B 



M - 



1 a 



D,, 



N 



A, 



(27) 

(28) 
(29) 



Please notice that after the steps following Eq. ([22|) the contributions of the non-diagonal ii^-terms are now contained 
in the doping dependent coupling constant in front of the last term of Eq. (^). At this stage we have also introduced 
a mass M for the bosonic variables which is arbitrary since the corresponding mass term M 2 ZZ gives only a constant 
in the path integral due to the constraint for the bosons fl2| . Resubstituting the functional derivative with respect to 
the source field in the fermionic part by the gauge field A M and expanding the remaining fermionic part and the 
bosonic one in powers of 1/ \N leads to 



(30) 
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In the limit of large N only the first two terms will survive. The first term is given by 



= ia(p = 0, co = 0) 



g/ff 



d 3 q 



1 



q 2 + v 2 + A/ 2 J ' 



(31) 



where v = 2irm/ (3 are bosonic Matsubara frequencies. We want to stress here that without taking into account the 
non-diagonal terms oc there would appear no doping dependent parameter as /* in Eq. (|3l]) since the diagonal 
fermionic contributions to are identical to zero due to the vanishing trace of the Pauli matrix a z . The second 
important term in the large A-expansion is given by the polarization insertions of the fermions and bosons for the 
gauge field A M as well as for the constraint field a: 



S{2) = \ J d3p {-^(P'^") ( n /L(P' w «) + 11 »u(P,u n ))\v(-p,-Un) - a(p,u n ) F(p,u n ) a(-p,-u) n )} 
where the polarization tensors are defined as follows: 



Q 2 + v m + M2 
2 



(q + p) 2 + (v m + u n ) 2 + M 2 
(2p + q) r (2p + q) s 



(32) 



(33) 



F(p,uj n 



V'i \ 

-{- 

2 \2m 

d 3 q 



A 77 



1 



i iv ™ - ^P 2 ][i(Vm + U) n ) - ^(p + q) 2 ] 



2m' 



<W / <f g 



1 



k - ^P 2 ][i{v m + w„) - 5^(p + q) 2 ] 
1 



q 2 + + M 2 (q + p) 2 + (u m + uj n ) 2 + M 2 



(34) 
(35) 



As a consequence of the local U(l) gauge-symmetry the polarization tensor II M „ = U^ v + U^ v has to be transversal 
p^Hfiv = and this property turns out to be very useful in evaluating the polarization tensors by means of a tensorial 
decomposition. 



A. From the renormalized classical to the quantum disordered phase by doping 



Since in the limit N — > oo the expression exp{— ^/NS^} gives a strongly fluctuating factor in the integrand of 
the generating functional (^) only the contribution = will survive. Evaluating the integral in the expression 
( pl| ) this saddle point equation establishes a definite relationship between the mass M of the Z-bosons, the coupling 
constant g, the temperature T, and the doping 6: 



M = — arcsinh 





■ exp 



-4tt/3 



■I r 



(36) 



with g c = 87r/ A denoting the critical coupling constant which depends on the ultraviolet cut-off A. Eq. (|36[) shows that 
the previously arbitrary mass M is now fixed by the saddle point condition and is therefore dynamically generated 
|^2| . By introducing the renormalized spin stiffness p = 2/g — 2/g c at zero doping, where the system is in the Neel 
ordered phase (at T = 0) we get the following expression for the mass M which is essentially the inverse of the 
spin-spin correlation length £ 



M 



£ 1 = —arcsinh 

P 



iexp 



{-2tt/?(p - 4 7 (1 + 2 7 m)S\ } 



(37) 



According to Chakravarty, Halpcrin and Nelson (CHN) |T3fl we have to distinguish three different regimes depending 
on the value of the argument y — ^ (p — 4 7 (1 + 2 7 m)d)~m the exponential of Eq. (p7|). For y ^> 1 the system is in 
the renormalized classical (RC) regime where the spin-spin correlation length £ diverges exponentially as T tends to 
zero, indicating a ground state with long-range antiferromagnetic Neel order at T = (/3 = oo) 



G 



£ = /3exp 



7r/3(p-4 7 (l + 27m)<s) . (38) 



Contrary for y -C — 1 the correlation length tends to a constant value at T = 

1 



4tt[47(1 + 27 m)<5 - p] ' 



(39) 



implying that the system is in a QD spin-liquid state. From the above discussion it is clear that there must be a 
phase transition due to doping at T — where the system goes over from the Neel state to the QD spin liquid state. 
The corresponding critical doping S c , where this quantum phase transition occurs, is defined by: 

(40) 



47(1 + 2 7 to) 



Thus the critical doping 5 C is completely determined by the experimental value of the spin-stiffness p at zero doping 
and the microscopic parameters m and 7. For an estimate of the critical doping S c we use a parameter set |lj| of 
the underlying three-band Hubbard model, which lies in the strong coupling regime, where the spin-fermion model 
applies. Also this parameter set is very realistic as was shown by quantum Monte Carlo simulations || through a 
comparison of numerical results and experimental data for the cuprates. Furthermore the assumption of isotropic 
masses and coupling constants in the theory made in the introduction is justified using this parameter set. Using the 
expressions for the mass to and the coupling constant 7 given in Ref. || we get m = 3/2 and 7 = 1/12 |^5[. With 
the experimental value p — 1.7 • 10~ 2 [|l6| we explicitly get S c ~ 4% plj . 

Finally for \y\ -C 1 the system is in the quantum critical (QC) regime which is completely determined by the critical 
point (S c ,T = 0) 0. 



B. Propagator of the gauge field 



After we have discussed the phase diagram of the doped AF we examine here the second term of the large- TV 
expansion, which consists of the polarization insertions of the fermions and bosons for the gauge field A M and for 
the constraint field a. As already stated, the polarization tensor ( |32"| ) of the gauge field A M has to be transverse and 
therefore the following tensorial decomposition is useful [|L7| : 



n^i» = [rf (p) + rf (p)] + b^( p ) [if r» + if (p)] , 



(41) 



with p = (p, w n ) denoting the wave- vector p together with the Matsubara frequency uj n . Here we have defined the 
two transverse tensors (p 2 = p 2 + ui 2 ) 



_ P^Po \ P^_ 
p 2 J p 2 



<W - 



PoPi- 
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(42) 



with the following properties: 



A,,,, -f- B n 



Suu - ^ ; A UV B^ = 



(43) 



The scalar functions T^' B (p) and r^' B (p) are evaluated in Appendix |A] and are given by: 



r» : =rf( P ) = rf(p) = 
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f p) 2 - i [v m + LO n ) 



(44) 
(45) 



(46) 
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After an analytical continuation to real frequencies iuj n — > u> + i0 + and evaluating the various integrals at T = we 
get the following expressions for the scalar Y functions (The details of this calculation are delegated to Appendix |^) : 

(47) 

(48) 
(49) 

where x = Pf is the fermionic density and vp the Fermi velocity. Within the Coulomb gauge, V ■ A = 0, we 

get for the propagator of the gauge field A M 




D^ = ^(r B (p,u J )+T[(p^))- 1 , (50) 
Doi = - (5!) 

B % = (% - P -§f) (r B (p,c) + rf( P ,c))- 1 . (52) 



Since we are interested on the long-wavelength low-frequency physics of the gauge theory it is appropriate to study 
the limit p — > 0, u> — ► 0. Because of the step functions in the fermionic contributions we have to distinguish two 
different regions in the ui-p plane. 



In the region, where x > 1, we get 



A > _ uj 2 -M 

TT — ~J2 i n . n„M ' V 53 ) 



P 2 -1-UJ 2 + BEM. 
1 24tt ~ m 

U V - [ °V „2 ) _J_..2_ PfM 



and for x < 1 we get: 



D rr = ^ p2 + ^ • (55) 

(56) 



D ij (p,«)= * 



P 2 ) l p 2_ ipF M 



PVF 



The propagator of the constraint field a is completely determined by the bosonic polarization and is therefore inde- 
pendent of the value of x 

D a = Fg 1 = 8ttM . (57) 

As we see from Eqs. (|5^) - ( |57j ) the time component of the gauge field propagator as well as the propagator of the 
Lagrange-multiplicr field a are massive, since they are Debeye screened, and thus mediate only short-range density- 

2 

density interactions between the fermions and bosons. Here we want to remark that the pole in the prefactor in 
Eq. ( p3| ) is not physical, since it is a consequence of the choosen Coulomb gauge. In fact one can explicitly see in the 
axial or in the Lorentz gauge, that the singularity due to that prefactor in Eq. ([33]) is absent. Physical poles of the 
gauge field propagator are only given by the zeros of the gauge independent combinations of the fermionic and bosonic 
scalar functions (T B (p, oj) + Tf(p, uj)) and (T B (p, u>) + r^(p, oj)) appearing in the corresponding expressions for the 
propagator. 

In the region |pjy^| > 1 the spatial component of the gauge field propagator has a pole at a finite frequency 

lu 2 = 247r PF J^ 1 with a form reminiscent of the plasmon pole in an electron gas. Furthermore, in the opposite limit, 
Ij^jjrl < 1, the gauge field propagator turns out to be massless and describes an overdamped mode with a dispersion 

law oj oc ip 3 , whose physical origin can be traced back to the Landau damping oc iuj/p due to fermions near the Fermi 
surface. As a consequence of this massless mode the fermions and bosons experience an effective long-range interaction 
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by exchanging transverse gauge bosons. The form of the propagator in Eq.(J5q) is essentially the same as found by 
Nagaosa and Lee [[l8| and Reizer Jl9|, with a central difference given by the fact that the mass M of the magnetic 
excitations determines the strength of the propagator. Hence as a special feature of the present theory a new energy 
scale, the spin-gap A cx M measured in neutron scattering experiments [^p| , enters in a natural way and appears as 
the strength of the interaction. Therefore, we expect that the long-range current-current interaction mediated by the 
transverse gauge field becomes only effective in the QD phase, where the spin-gap A is large compared to the thermal 
energy ksT. This implies that the tendency to form bound states only occurs in the disordered magnetic state, such 
that the phenomena of charge-spin separation and pairing in this theory should be restricted to that region. 



III. DISCUSSION AND SUMMARY 



As was noticed first by Holstein et al. 21 and later by Reizer figfl , the coupling of fermions to a transverse vector 
gauge field leads to deviations from the Fermi-liquid picture. However, in the case of a conventional metal interacting 
with the electromagnetic field, the effects of the exchange of transverse gauge bosons are in practice unobservable since 
they are suppressed by the ratio of the Fermi-velocity and the speed of light (— ) ■ In contrast to those systems, in 
the above gauge theory for the doped AF the coupling constant is of order unity and therefore no such small parameter 
enters, a feature common to gauge theories developed recently for strongly correlated metals |l^,|2^,^]. Calculating 
the self-energy for the fermions at k p in lowest order of perturbation theory one gets an unusual frequency dependence 

E oc uji , which implies that the quasi-particle weight Z — 1 1 — ^jRe E w ^o| vanishes at the Fermi surface signaling 
a break down of Fermi- liquid picture fig] , and hence, a perturbative expansion in the coupling constant is in this case 
not valid pq] and one has to implement some non-perturbative approaches. In spite of the conflicting claims in the 
recent literature concerning the one-particle properties of gauge theories for strongly correlated metals ]24[| , we would 
like to point out that the results obtained there are not directly applicable to our case, since the physically relevant 
one-particle properties appear only after a rotation back to the global refence frame. This amounts in the lowest 
order of the large TV-expansion to a convolution of the free fermionic and bosonic propagators and the corrections 
due to the gauge fields are obtained at 0(1/N). Results pertinent to this point will be presented elsewhere |p5| . An 
even stronger difference with the general discussion of gauge theories in this context should be expected due to the 
possibility of formation of bound states, that we discuss in the following. 

The singular mode of the gauge field mediates a long-range current-current interaction between the Z-bosons and 
the fermions, which in the static limit reduces to a logarithmic potential. Moreover since the strength of the singular 
gauge field propagator is proportional to the spin-gap A we expect the formation of bound states only in the QD 
phase, where the spin-gap is large compared to the thermal energy ksT . Hence particles with opposite charge will 
bind and the low-energy sector of the physical spectrum of the theory will contain only states with zero charge with 
respect to the gauge field A. Since the charge of the fields is given by their spin-projection, the bound states correspond 
to spin-singlets. According to table ( |l3|) the following three bound states are possible: 

(a) Bound states between two bosons: 

Z-Z (58) 

These bound states correspond to spin-waves excitations (remember n = ZcrZ) at the antifcrromagnetic wave 
vector Q — (tt, it) with a gap A cx M in the spectrum. 

(b) Bound states between a fermion and a boson: 

Z-p (59) 

These bound states have an electric charge +e but have no spin. Thus they describe spinless charge-excitation 
with fermionic character. Thus, this scenario gives an alternative way to charge-spin separation, where the bare 
excitations are just spin-i fermions but the renormalized ones are spinless. 
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(c) Bound states between two fermions 



Pi - .Pt ( 60 ) 



Such a bound state describes the pairing of two fermions in a spin-singlet with charge 2e. Hence, charge-spin 
separation and pairing are intimately connected in our case and result from the same interaction. 
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A similar description for holes in a spin-liquid state was derived by Wen [g6| starting with a one-band Hubbard model. 
However, in that case charge and spin degrees of freedom are separated from the very beginning (the bare fermions 
do not carry spin) and hence the phenomenon of charge-spin separation does not follow from the formation of bound 
states as in our theory. Quite on the contrary, the bound states between fermions and bosons in that theory carry 
charge e and spin, implying that charge-spin separation only occurs when there is no long-range interaction present, 
i. e. when either the gauge field acquires a mass due to a Chcrn-Simons term or abnormal screening at a high doping 
level takes place p7| . 

For the bound states (|5|)-(^0|) we get the following picture. At a doping level which is greater than the critical 
one, there is a crossover from the QC phase at high temperatures to the QD phase at low T. Below the crossover 
temperature T* the bound states will form and this should be reflected by various physical quantities. 

The fact that the formation of the bound states brings the system towards a singlet state, may explain the anomalous 
data in nuclear magnetic relaxation experiments pq| , where a strong reduction of both the Knight-shift and the 
relaxation-rate on oxygen and copper is observed as a function of temperature, well above the superconducting 
transition, for underdoped samples of YBCO. Closely connected with this is the observation of a spin-gap in neutron 
scattering experiments on YBCO [^o| which corresponds to the mass of the Z-bosons. Moreover, when the bound 
states between a fermion and a boson form, charge-spin separation takes place and the scattering of the charge- 
carriers by the gauge fields vanishes since the new bound state has zero charge and therefore does not couple to 



the gauge field. Such a scenario is consistent with recent resistivity measurements in underdoped YBCO 123 |E9 30 
and 124 f3lf| , where a reduction of the resistivity with respect to the linear temperature dependence is observed, at a 
temperature T* where the nuclear relaxation rate (1/TiT) on Cu shows a maximum that is commonly associated with 
the spin-gap p8|. Furthermore in the resistivity perpendicular to the CuO-planes a crossover from a metallic-like to 
a semiconductor-like behavior is observed at the same characteristic temperature T* p0| , which in our theory follows 
from the fact that one has to break the bound states in order to have charge transport perpendicular to the planes. 
This is also consistent with a suppression of the optical conductivity parallel to the c-axis that was reported recently 
for underdoped samples of 123 |32] indicating the opening of a pseudo gap along the c-axis. 

In conclusion we have shown that within the long-wavelength, low-energy sector of the spin-fermion model a 
quantitative correct description for the doping induced quantum phase transition appearing in the HTSC materials 
can be obtained. Furthermore, in the gauge field theory of the doped AF in the rotating reference frame, charge-spin 
separation and pairing follow in a natural way from the formation of bound state in the QD phase (above the critical 
doping and at low temperatures) due to the long-range force mediated by the singular mode of the transverse gauge 
field, such that an intimate connection between these phenomena and the spin-gap results. Finally, we have also 
discussed the connection of these bound states with the anomalous features observed in the low-doping regime of the 
cuprates. 



We acknowledge support by the Deutsche Forschungsgemeinschaft under Project No. Mu 820/5-2. 

APPENDIX A: TENSORIAL DECOMPOSITION 

(a) Evaluation of rf 

Due to the properties ( fi"3"| ) of the two transversal tensors A^ v and we get (p 2 — p 2 + u 2 ) 

('*> - P ~f) £ (** - P ~f) KM = - rf (p) , (Ai) 

where Hoo(p) follows from the expression ( |34] ) of the fermionic contribution to the polarization tensor: 

n °°(p) = \ E / 7^7^ TT- " - 2 T =: H(p) • (A2) 



(b) Evaluation of rf 

From (f43f) further follows: 



*» (*« ~ P -f) **> K = - P ~f) K = n£ - v -fK = rf , (A3) 
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with 



m (5 2-~> 



d 2 q 



1 

2^ 



d 2 q 



(2q + p) ■ (2q + p) 



2i(9+P) 2 -*("m+Wn) 



(A4) 



and 



1 2lv- 



m 



d 2 ? 



( 27r ) 2 2b? 2 



i \ 2 1 

2m) 



-:^f(p) 



d 2 q 



(2q + p)-p (2q + p)-p 



1 

2^i 



(A5) 



APPENDIX B: FERMIONIC AND BOSONIC INTEGRALS 



In this appendix we calculate the bosonic and fermionic integrals appearing in the scalar functions T , rf and T 



of the polarization tensor for the gauge field. 

First we evaluate the bosonic integral F(p,uj n ) defined by: 

d 2 q 



1 



r) 2 [q 2 + v 2 n + M 2 ] [(q + p) 2 + (y m + oj n ) 2 + M 2 } 
Using a Feynman parametrisation we obtain: 



M 2 + q 2 + v 2 + p 2 x(l -x)+ uj 2 x{1 - x) 

where v = v + w(l — x). Integration over momentum leads to: 

^ If 1 1 

F (P^n) = -g ]T ^ J dX M 2 + (p 2 +U 2 )x(l-x)+V 2 ■ 



(Bl) 



(B2) 



(B3) 



In order to evaluate the frequency sum we make a partial fraction decomposition of the integrand. (Here we have to 
introduce a factor exp(ivr]) to ensure convergence of the series): 



F(p tUn )=^J o 



-^2^/M 2 + (p 2 +Lo 2 n )x(l-x) 
1 



exp(ivrj) 



iv + [^JM 2 + (p 2 + ujI)x{1 -x) + iu n {l - xfj iv - [^M 2 + (p 2 + ul)x{l - x) - iu> n (l - i)) 
1 dx 1 



1 + 



ii 4tt 2y/M 2 + (p 2 + u:l)x{l - x) 
1 



(B4) 



exp 



Py/M 2 + (p2 + u 2 ^x{\ -x)- ipu> n {l - x) 



+ 



— 1 exp 



Py/M 2 + (p2 + uj2)x(1 -x)+ i(3u> n (l - x) 



1 



For T — > the temperature dependent part of the integral gives only an exponentially small correction and can 
therefore be neglected. After an analytical continuation to reel frequencies ioj n — > ui + irj we get: 
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F(p,u>) 



T^O 1 



47T _ w 2 



arcsm 



P 2 — LO 2 



AM 2 + {p 2 - to 2 ) 



(B5) 



In the following we calculate the fermionic integrals. The first one is gives and equals the density of fermions at the 
fermi surface: 



1 



d 2 q _ _ 



= ~PF ■ 



For the calculation of the remaining fermionic integrals we make use of the formal identity 



1 



V- T «7r S(lo) 



LO ± IT] LO 

with V denoting the principal value. In this way we obtain 

fiZsJ (2^)2 [_i_ (7 2_ Wm ] \j_ i[q+p) 2 _ Wm _ [uJ + l1l) 

d 2 q n(q - p/2) - n(q + p/2) r ^ 2 



(B6) 



(B7) 



- [#1 

J (27T) 



[n(q - p/2) - n(q + p/2)} 5 (w - (q ■ p) /m) . (B8) 



r) 2 LO — (q ■ p) /m 

Furthermore, we go over to zero temperature, where the difference between the two Fermi distributions is given by 

n(q - p/2) - n(q + p/2) = p ■ Vn(q) T =° p ■ q/q 5{\q\ - k F ) . (B9) 
Using this expression we get 

Ii = -V(±f2 fdB rdqS(q-k F ) ; PCOS ^, 

2ir J J Q to - (qpcos9)/m 

1 f 71 f°° 

+ i — d9 dq qp cos 9 S(q — kp)5 (lo — (qp cos 9) /m) 
27r Jo Jo 

=L= + i(^)m J dzS(z-x) 



2m — 
2tt 



dz 

1 Z — Xy/1 



m 
2^ 



xO(x 2 -l) xQ(l-x 2 ) 
1 1==^- +i- 





LO 


) ; x = 






pv F 



\/x 2 — 1 \J\ — X 2 

The remaining fermionic integrals are calculated in the same way and we just give the results here 

(2q + p) ■ (2q + p) 



(BIO) 



Iv f dq 



IV, 



»] 2k + P) 2 - iv m - (w + if?) 



and 



4fc|/ 2 F 



1 y 



(Bll) 



d 2 q 



(2q + p)-p (2q + p)-p 



(2n) 2 r 

L 2m " 



1 „2 



If, 



»] sir (9 + p) 2 - - (w + if?) 



-4fcV ( ^ + - 



1 x 3 e(a; 2 -l) x 3 6(l-x 2 ) 



(B12) 
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